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We derive the scaling structure of the Kosterlitz-Thouless-Berezinskii (KTB) transition temper- 
ature of a homogeneous Bose gas in two dimensions within diagrammatic perturbation theory. 
Approaching the system from above the transition, we calculate the critical temperature, Tkt, and 
show how the superfluid mass density emerges from Josephson's relation as an interplay between the 
condensate density in a finite size system, and the infrared structure of the single particle Green's 
function. We then discuss the trapped two-dimensional Bose gas, where the interaction changes the 
transition qualitatively from Bose-Einstein in an ideal gas to a KTB transition in the thermody- 
namic limit. We show that the transition temperature lies below the ideal Bose-Einstein transition 
temperature, and calculate the first correction in terms of the interparticle interactions. The jump 
of the total superfluid mass at the transition is suppressed in a trapped system. 
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I. INTRODUCTION 

The ability to produce a two-dimensional atomic gas 
trapped in an optical potential has stimulated interest in 
the Kosterlitz-Thouless-Berezinskii (KTB) transition in 
such a system 0, J, 3, 4]. Recently, interesting phase 
patterns have been measured in a quasi two-dimensional 
situation j^. A homogeneous Bose gas in two dimen- 
sions undergoes Bose-Einstein condensation only at zero 
temperature, since long wavelength phase fluctuactions 
destroy long range order ; nonetheless interparticle in- 
teractions drive a phase transition to a superfluid state at 
finite temperature, as first pointed out by Berezinskii 
and by Kosterlitz and Thoulcss 8] . The phase transtion 
is characterized by an algebraic decay of the off-diagonal 
one-body density matrix below the transition tempera- 
ture, Tkt- Furthermore, the superfluid mass density, ps, 
jumps with falling temperature, T, from just above 
Tkt to a universal value ps = 2m'^TKT / just below 
where m is the atomic mass. We use units h — ks — 1- 

In a non-interacting homogeneous Bose gas in two di- 
mensions, the density is given in terms of the chemical 
potential p by 



log[l 



(1) 



where f3 = 1/T, and A = (27r/mT)i/2 jg ^t^^ thermal 
wavelength. As p approaches zero, the phase-space den- 
sity grows arbitrarily, implying that the ideal Bose gas 
does not undergo Bose condensation at finite tempera- 
ture. On the other hand, a trapped non-interacting gas 
undergoes a Bose-Einstein condensation at a finite tem- 
perature |lO|. In the semi-classical limit, the total num- 
ber of particles is given by N{fip) = g2{—Pp){T/Lj)^, 
where 52(2;) = X^J^i e~^^/n'^ ~ 7r^/6 — x{l — logs), for 



a; — > 0, and uj is the trapping frequency; thus the con- 
densation temperature is Tbec — V&Nlo/t:. 

The question of whether an interacting trapped Bose 
gas undergoes an ordinary Bose-Einstein or a KTB has 
not completely been settled j^J. In this paper we show 
that in the thermodynamic limit - in which the total 
particle number N goes to c», and the trap frequency, 
Lu, goes 0, with Nlu'^ constant - interactions at the mean 
field level destroy ordinary Bose condensation; instead, 
the system undergoes a KTB transition at a temperature 
slightly below the ideal condensation temperature, Tbec- 

We approach the phase transition of the homogeneous 
two-dimensional Bose gas by carrying out a scaling anal- 
ysis (Sec. II and III) similar to the one we used earlier to 
discuss the phase transition in a dilute three-dimensional 
homogeneous Bose gas El El El- The phase below TifT 
is characterized by an algebraic decay in space of the sin- 
gle particle Green's function, G. As in three dimensions, 
the transition point occurs when the single particle spec- 
trum becomes gapless, G'~^(0, 0) — 0, where G{k,z^) is 
the Fourier transform of G in space and time. We first 
discuss, in Sec. II, the scaling structure of G as Tkt 
is approached from above, and rederive the relation be- 
tween the temperature and density at the transition: 



27r 



Tkt 



n mlog{G/a)' 



(2) 



where C is a constant 1151 . ~ (380 ± 3)/27r according to 
numerical simulation [la |. and we take the interparticle 
coupling constant in two dimensions to have the form 
g = 2TTa/m, where a is a dimensionless parameter 21]. 

We then turn, in Sec. Ill, to the structure below Tkt, 
working in terms of the Green's function in momentum 
space. Below Tkt we expect G{r — r') ~ l/|r — r'|'' for 
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an infinite system and |r — r'| —^ oo, where 77 depends 
on T (and equals 1/4 at T^t @); however, G{k) is not 
well-defined in an infinite system. We therefore analyze 
the structure starting with a finite size system, of charac- 
teristic dimension L; such a system has a non-zero con- 
densate density, no ^ G{r ^ L) ~ l/i^', which vanishes 
in the thermodynamic limit, L 00. Nonetheless, we 
show using Josephson's relation between the condensate 
density, the superfluid mass density, ps, and the infrared 
behavior of G{k) ^M,E^, 



(3) 



that ps is finite in the thermodynamic limit. Our analysis 
of the scaling structure is valid for all coupling strengths, 
and thus, going beyond perturbation theory, is an exten- 
sion of the work of Popov ^] . 

Then, in Sec. IV, we apply our results to a trapped 
gas within the local density approximation. We show 
that the KTB transition temperature, for a weakly inter- 
acting system, lies below the Bose-Einstein condensation 
temperature of the ideal trapped gas by terms of order 
Qflog^a. Furthermore, the jump of the total superfluid 
mass at the transition is a, and is thus highly sup- 
pressed compared with that in a homogeneous system. 



II. SCALING STRUCTURE ABOVE THE 
TRANSITION 

In this section we derive the KTB transition of a two- 
dimensional weakly interactiong homogeneous Bose gas 
by studying, as in 0, the scahng structure of G just 
above the transition. For wavevector k and complex fre- 
quency z, G{k, z) is given in terms of the self-energy, 
S(/c,z), by 

G-\k,z)^z + p-^-^{k,z), (4) 

where p is the chemical potential. The transition is 
defined, for fixed density, by the temperature at which 
the single particle spectrum becomes gapless, and conse- 
quently off-diagonal single particle density matrix decays 
algebraically; at the transition point, G'~^(0,0) = 0, or 
p ~ S(0, 0). In terms of G the density is 



(pk 



:G{k, Zi^). 



(5) 



where the Zi, = 27rwT are the Matsubara frequencies 
(i/ = 0,±l,±2,...). 

The mean-field contribution, = 2gn, to the self 

energy, is independent of frequency and momentum, and 
can be absorbed in a shift of the chemical potential. We 
introduce the mean-field coherence length, ^, by 



1 



^mf ^ P\ 



(6) 



at the transition, 

l/2mC' = E„/~S(0,0). (7) 
In terms of C, the mean-field Green's function is, 

2toC^ 



Gmf{k,Zi,) 



(fcC)2 + l-8i7r2(C/A)2jy' 



(8) 



In the zero Matsubara frequency sector, the self-energy 
has the scaling structure: 



(9) 



where ctq is dimensionless and ~ as a — > 0. In contrast 
to the case of three-dimensions, self-energy diagrams be- 
yond mean field are ultraviolet convergent even in the 
zero Matsubara frequency sector. 

Using Eqs. ijHl and © we may write the transition 
condition, ([TJ), as 



A2 



<To ( 0, ) + 1 = 0; 



(10) 



at the phase transition the parameter a(^^ /X^ = J ap- 
proaches a finite value(fixed point), J*, determinined 
by this equation. The contribution of non-zero Mat- 
subara frequencies in the denominator of Eq. is ~ 
{a('^ / X^) / a, and thus in the limit of small a, the con- 
tributions of non-zero Matsubara frequencies to the self- 
energy are of relative order a and higher, and can be 
neglected. 

To calculate the critical density at given temperature 
we use the mean-field density as reference: 

nc{a,T) = nmf(^ -l/2mC^) 
d^k 



(2^)2 



[G{k,z^) - Gnii{k,z^)] 



(11) 



In leading order, we neglect non-zero Matsubara frequen- 
cies in the summation on the right, and derive 



X^Uc = log + 2 I kdk 



1 



fc2+C-2(l+C7o) 



k^ + C 



Oia) 



(12) 



Using Eq. ^ we see that the integral in H12|l is a constant 
of order unity, independent of a. Thus we arrive at the 
critical density, 



\''n, ^ log- +0{a), 
a 



(13) 



as in Eq. 

The true correlation length, ^, above Tkt is given by 



2m^ 



^ = S(0, 0) - M = S(0, 0) - S™/ + — 



2mC2 



(14) 
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To make contact with the theory of critical phenomena, 
we rewrite this equation as 



M{J) 



1 
J 



2toA2 



(I](0,0)-E„/), (15) 



where from the previous discussion we know that the 
right side is a function only of J = a^^/A^. At the tran- 
sition, J goes to the fixed point, J*, given by M{J*) = 0. 
Furthermore, approaching the transition from above, we 
see from Eq. ^ that A J = J - J* cx (T - TKT)/a. 
The dependence of M{J) on A J near the fixed point de- 
termines the critical index v, of the correlation length, 
f - |r - TktI'", in particular that M(J) - jAJp"^ - 
|T — Tkt\'^^ , in the neighborhood of the fixed point. 



III. SCALING STRUCTURE BELOW THE 
TRANSITION 

At the transition temperature, T^k, the single particle 
Green's function decays algebraically. Below, the scaling 
structure is most readily analyzed in momentum space, 
as above; however this approach is made difficult by the 
fact that below Ttk the single particle Green's functions 
continues to decay algebraically in real space sufficiently 
slowly that its Fourier transform is not absolutely con- 
vergent. In order to avoid this problem, and to use the 
same approach as above, we adopt the strategy of work- 
ing in a finite size system, of characteristic dimension L, 
in which the condensate density, no, is non-zero. At the 
very end, we take the thermodynamic limit, L — > cxi, at 
fixed density, in which case no also goes to zero, since 
long-range order is prohibited in two dimensions. 

In a finite size system the condensate density is given 
by no = n — fi, where 



G{k,n) 



(16) 



and fco ^ 7r/i; in leading order, G here can be taken to 
be the infinite size Green's function. Since for T < Tkt, 



G{k -> 0,0) = ~2mKC/k 



2-r, 



(17) 



where the constants K and ?7 (0 < 77 < 2) in general 
depend on T, we find to leading order. 



(18) 



The resuh ^ is given for a circular box; in general it 
is modified by a numerical factor close to unity, weakly 
dependent on the geometry. 

Below Tkt the system is superfluid even though the 
condensate density is not extensive and vanishes in the 
thermodynamic limit. However, in two dimensions, an 
algebraically decaying correlation function is sufficient 
to yield a non- vanishing superfluid density, ps, in the 



thermod yna mic limit , as can be seen from the Josephson 
relation \vA llsL l25l | between pg and no in a finite 
system: 



lim 



nQin 



k~.k„BG{k,Qy 
Thus from Eqs. (Ejl and and find 



Ps 



2771] 



(19) 



(20) 



Inverting the relation we obtain 77 = m?T /2t:Ps] thus this 
prediction of the spin- wave approximation for the critical 
exponent is prediction below Tkt- Since ps < mn, we 
see that 



> 



1 



(21) 



Furthermore, since t] < 2, the superfluid mass density 
can never vanish at finite temperature below Tkt, and 
thus it must be discontinuous across the phase transi- 
tion. At Tkt, the value of the universal jump of the su- 
perfluid density predicted by Nelson and Kosterlitz j^, 
p^ = ATKT'm'^/^TT, leads to = 1/4. 

Let us now discuss the detailed structure of the Green's 
function below the transition. We basically follow the 
scaling approach used in j23j | in three dimensions. Our 
strategy is to expand the self-energies formally in powers 
of a, no, and fcg. In the infinite size system, kg — s- 0, the 
self-energies diverge as no 0; the point, fco — + 0, no — > 
is singular. However Joscphson's relation constrains the 
limit no — > and fco — > in terms of ps- 

The particle density, n, in the condensed phase 
is a function of a, no, and T, and has the form, 
n(a, no, ko,T) = no + n(a, no, fco, T), where n(a, no, T) is 
the density of non-condensed particles (with momentum 
k > ko). At the transition temperature, n(a,0,Tc) = Uc- 
We calculate n(a, no,T) in terms of the matrix Green's 
function. 



g{rt,r't') = -i{{T{^{rt)^\r't'))) 
-{^\r't')){^{rt))), 



(22) 



where the two component field operator is ^'(rt) — 
{^p{rt),tp'^ {rt)^ The Fourier components of Q^^ have 
the form, 



G ^{k,Zn) 



_ ( Zn + p- ek-T'ii 



-'21 



-Zn + p — £k 



^22 



(23) 



where Sk — fc^/2m, and the Yiij{k,Zn) are the corre- 
sponding self-energies. The chemical potential, p, de- 
pends on no and fco, and is specified by the Hugenholtz- 
Pines relation Ejl, 



// = En(0,0)- £12(0,0). 



(24) 
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The lowest order mean field self-energies, En 



■"11 



2g{nQ + h), S12 = = 5*^0 are independent of mo- 

menta and Matsubara frequency, and, as above Tkt, we 
absorb them in a mean field coherence length, ^, 



mf 



(Eii(0,0)- 
-1/2toC^. 



Ei2(0,0))-(I]r/-E 



12 ) 



(25) 



Since the propagators remain formally infrared con- 
vergent we can derive the scaling structure of the self- 
energies by power-counting. As above the transition, we 
may neglect non-zero Matsubara contributions to lead- 
ing order. As opposed to the three-dimensional case, 
there are no formal ultraviolet divergencies in the ex- 
pansion beyond those in mean field, and therefore need 
for renormalization. The expansion of the self-energies 
beyond mean field starts at order a^(^^/A^; furthermore 
S12 is formally at least of order uq. Diagrams of order 
g'^ with K > 3 in the formal expansion contain vertices 
with two Green's functions entering; similar to the struc- 
ture at Tc, they involve the dimensionless combinations 
ctC^ jy? = P and uqX^. The latter part originates from 
the dependence of G"^-^ on 2mY:I^2 ~ oluq. Any dia- 
gram with an explicit power, p, of no can be generated 
from a corresponding diagram of power p — 1 in which a 
line is replaced by at each of its ends. Thus each 
power of no involves one fewer two-momentum loop to 
be integrated over. The explicit no dependences enter in 
two ways. Terms involving G\\ — G12 lead to the com- 
bination P^uqX^, which vanishes as no — > 0. On the 
other hand terms involving the combination Gn + G12, 
which in mean field diverges as no — s- in the infrared 
limit, lead to divergences which are cutoff by fco and 
thus produce an additional (koC)'^ / {hqX^ P) dependence. 
In the limit no — *■ 0, fco — > 0, only the dependence on 
[koC)'^ / uqX^ = Q survives. 

Then with all momenta fc scaled by 1 /^, we find the fol- 
lowing scaling structure for the self-energies in this limit. 



2to 



(e,, (fc) - Ei;;^ (0)) - a,, ikC,P,Q), (26) 



where the aij are dimensionless functions of dimension- 
less variables, and we neglect terms proportional to pos- 
itive powers of no and fco. In particular, for fc — > 0, 



2m 



((En(0,0)-Ei2(0,0))-(E5'1^E1'^0) 

= s(P,Q), (27) 



where s is a dimensionless function. Comparing with 
Eq. H25|l we see that s -I- 1 = 0, an equation that deter- 
mines P as a function of Q, so that 



c = ^/i(g). 



(28) 



where h is dimensionless. 



We now take the limit no — > and fco 0. Equa- 
tions (Uni imply that Gii(fco,0) = -m(^£{Q), and 
thus the Josephson relation implies that the superfluid 
mass density has the structure. 



(29) 



A2Q£(Q)' 



where £ is a dimensionless function; this equation relates 
no and fco to ps in the limiting process. In this limit P 
as well as Q are constants dependent only on Ps/T. Sim- 
ilarly the other dimensionless functions depend in this 
limit only on ps- These forms are valid over the tem- 
perature range from Tkt down to T = (with non-zero 
Matsubara frequency contributions neglected); this en- 
tire region is therefore critical. 

We can gain further insight into the structure be- 
low by reformulating it in terms of the correla- 
tion lengths, Ct and C^, that control the infrared behav- 
ior of the transverse and longitudinal Green's functions, 
Gt — Gn — G12, and Gl — Gn + G12. We regularize 
the infrared divergent structure below Tkt by assum- 
ing a finite condensate density, no, and finite correlation 
lengths. In this description, the low temperature phase 
of the KTB transition is characterized by the ratio of 
amplitude (L) and phase (T) fluctuations of the order 
parameter, even in the absence of long-range order in 
the thermodynamic limit in which no ^ 0. In the end 
we take the limit, no — > 0, Ct ^ 00, and Cl ^ 00. In 
this way, we do not have to introduce an explicit infrared 
cutoff, fco, as we had to above. We define 



En(0,0)-Ei2(0,0)-/^, 



and 



— ^ = En(0,0) + Ei2(0,0)-Ai, 
zJrnC,^ 

in terms of which Gt and Gl are given by, 

fc2 



(30) 



(31) 



G^r 



2m 



+ [En(fc,0)-Ei2(fc,0)] 



-[En(fc,0)-Ei2(0,0)H- 



1 



2mCj 



(32) 



and 



Gii(fc,0) = - — + [En(fc,0) + Ei2(fc,0)] 
* zm 



[En(fc,0)-KEi2(0,0)] 



2™C? 



(33) 



The equilibrium state of the system is specified below 
the transition by Ct ^ 00, as seen from the Hugenholtz- 
Pines relation (|24|l . as well as no 0. From Eqs. H3U|) 
and H31|l . we obtain the ratio 



|-) =1 + 4toCtSi2(0,0). 



(34) 
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The second term is of order UqC,^. From the previ- 
ous discussion, we deduce that in the hmit uqX^ 0, 
4mCT5]i2(0,0) 8TTano(^f{CT/CL), where / is a di- 
mensionless function. Equation (|34|) thus implies that 
the ratio Ct/Cl is a function only of no(^. We there- 
fore we arrive at a self-consistent solution in the limit 
TT-o ^ and Ct — > oo in which uqCt remains finite. The 
value of this parameter is determined in terms o f by 
Josephson's relation, which we write in the form [l9j 



Ps — lim mno I 1 + 

2mA[En(fc,0)-5]i2(fc,0)]|,^J.(35) 

For finite (t, the derivative of the self energies with re- 
spect to are finite in the limit A: — > 0, but must diverge 
as Ct- The right side is a function only of noC-r 
limit C,x ^ oo, and no — > 0, thus defining uqCt in terms 

of Ps. 



IV. THE TRANSITION IN A TRAPPED GAS 

We turn now to the behavior of a two-dimensional sys- 
tem trapped in an oscillator potential, of frequency to. 
We consider for simplicity only the thermodynamic limit 
TV ^ oo, w — > 0, with Noj'^ constant, where N is the 
total particle number. 

In the absence of interactions the system undergoes a 
Bose-Einstein condensation at the critical temperature 
Tbec = V6Nlj/tt. However, at the mean field level, 
interactions destroy simple Bose-Einstein condensation. 
To see how the critical temperature is reduced to zero, 
we begin with the density profile calculated in the local 
density approximation, 



1 



(27r)2 e/SCfcVani+KttlO-p) _ i 



^^log(l-e^(^-^-«) 



(36) 



where 14fi(r) = -|- 2gn{r). Since n(r) has neg- 

ative curvature at the origin, interactions tend to reduce 
the effective trapping frequency. Explicitly, expanding 
(|36|l to order about the origin, we find that the self- 
consistent trap frequency at the origin is given by. 



25 



9^»mf(?-) 
1 _ g/3(M-2g«(0)) 



4/3.gA-2 + 1 - e/3(M-23n(o)) ' 



(37) 



As ^ 2gn(0), the limit in which Bose-Einstein conden- 
sation occurs in mean field, we have 



^eff ^ (25^(0) - p) ; 



(38) 



the potential becomes arbitrarily flat as /i — > 2gn(0). 
In the thermodynamic limit, Nuj^ constant, the ef- 
fective oscillator length d^s = 1 / y/mujcs is ~ 
{NT/{2gn{0) - p))^^^. At the point of mean-field Bose 
condensation, 2gn{0) — p = l/2mC(0)^ is therefore ^ 
T/N; using the estimate Q), we see that the num- 
ber of particles in the center, n{0)Trd'^g, grows as 
{T/uj)y/a{({0)/X)\og{({0)/X). Now, an extensive pop- 
ulation of the condensate mode implies C(0)/A ^ N^^"^ 
within mean field. From our estimate of the number of 
(excited) particles in the center of the trap we obtain 
an upper bound for the critical temperature of Bose con- 
densation, ~ ujN^I"^ I (o^l'^ log A'^). In the thermodynamic 
limit of the trapped system, the critical temperature of 
Bose condensation goes to zero. 

Mean fields destroy Bose-Einstein condensation of an 
ideal gas. However, the system does undergo a KTB 
transition in the thermodynamic limit, at a temprature 
calculable, to leading order in a, in terms of the mean 
field density profile, Eq. (|36|l . The KTB transition occurs 
when the chemical potential reaches the critical chemical 
potential, 2(0,0), calculated for a homogeneous system 
of the same density and temperature as in the center of 
the trap, or equivalently, Eq. (|13|l is satisfied. As we see 
from Eqs. Hll|) and H12|) in the homogeneous case, the 
critical density is given to logarithmic accuracy by the 
mean field density evaluated at the critical C given by 
Eq. {Tj); the corrections arise from critical fluctuations, 
which are, however, important only at small distances 
< \gn{^) — or 



where ^mLuT 



< 



1 



(39) 



In the thermodynamic limit the critical region produces 
corrections to the total number, AN/N ~ log a, which 
can be neglected, and thus, to leading order in a, we can 
calculate the transition from density profile, Eq. (|3fi|l . 
Expanding n^nfif) to first order in a we have, 



ni^{{r) = no{r) 



1 



A2 



where 



no(r)^-^log(l-e'^(^— ^--^Z^) 



(40) 



(41) 



is the ideal gas density. Integrating over space, we find 
the first correction to the total number of particles. 



/ 2^no(0) Y 
\ mu) J 

2 



N ~ J (i^rnmf(r) = J d'^rno{r) — a 

T2 f27Tlog{-pp) 
= 92{-f3p)— - a 



(42) 



where p = 2gn{Q) — l/2mC^. At the Kosterlitz-Thouless 
transition, p ^ a. To within logarithmic corrections we 
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may replace no(0) by —A ^ log(a), and g2{—f3^) by (72(0). 
Thus at given T , 



NKTiT,a)^NBEciT)--a 



2tt log a 
rriLoX^ 



(43) 



where N-qeqIT) — {-k"^ / (!>){T / ujY is the relation between 
particle number and temperature of the ideal Bose gas 
transition in the trap. Equivalently 



Tkt — Tbec 



T, 



BEC 



-alog^a + ©(aloga). (44) 



The existence of the Kosterlitz-Thouless transition just 
below Tbec in a trap involves physics beyond mean-field. 

In a KTB transition in a homogeneous system, the 
superfluid mass density jumps with falling temperature 
discontinuously at Tkt from zero to 2m'^TKT/'^- In a 
trap, however, the transition first occurs in the center, 
and extends over a region of size Tc, Eq. H39|) . The total 
superfluid mass, M^, is therefore 



Ps{r)cPr ^ TTps{r = 0)r^ 



so that 



M 



1 

N 



A2 

e 



a 



(45) 



(46) 



where M = mN, and P = a^^/A^ is a constant (fixed 
point) at the central transition point. The jump in Ms 
is thus highly suppressed, and goes to zero in the limit 
of an ideal gas. 

A key indication of superfluid behavior below the tran- 
sition would be creation of a vortex at the center of 
the trap, where the system first becomes superfluid, 
e.g., by cooling a rotating system through Tkt- In or- 
der to create a vortex at Tkt it is necessary that the 
vortex core, of radius ~ C, fit within the critical re- 
gion, of size ~ Tc, at the transition. From Eq. H39|l . 
C/^c ~ mujC,'^, so that in the thermodynamic limit, 
C/t'c ~ P/a^fN. Thus for large creation of a vor- 
tex within the critical region becomes possible. The 
critical rotation frequency, fie, for creation of a vor- 
tex is of order (1/mr^) log(rc/C), and therefore Vlc/uj ^ 
(CAe) log(rc/C) ^ [l/aVN) \og{aVN). 

A further probe of the state of the system below Tkt 
would be determination of the density correlations, as 
have been recently measured in Bose gases trapped in op- 
tical lattices . These correlations will depend strongly 
on the amplitude fluctuations, described by C,t Wh - 
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